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Abstract 

The distortion of the regular motion in a quantum system by its coupling to the continuum 
of decay channels is investigated. The regular motion is described by means of a Poissonian 
ensemble. We focus on the case of only few channels K < 10. The coupling to the continuum 
induces two main effects, due to which the distorted system differs from a chaotic system 
(described by a Gaussian ensemble): 1. The width distribution for large coupling becomes 
broader than the corresponding Xk distribution in the GOE case. 2. Due to the coupling to 
the continuum, correlations are induced not only between the positions of the resonances but 
also between positions and widths. These correlations remain even in the strong coupling 
limit. In order to explain these results, an asymptotic expression for the width distribution 
is derived for the one channel case. It relates the width of a trapped resonance state to the 
distance between its two neighboring levels. 



1 Introduction 



Resonances are states of an open quantum mechanical system, in which the internal dynam- 
ics of the underlying closed system is disturbed by the coupling to the decay channels. For 
small coupling strength the widths of the resonances increase with growing coupling while their 
positions remain almost unchanged [l], ^, (3|. Therefore the degree of overlap of neighboring reso- 
nances increases. This gives rise to interferences and the internal dynamics may suffer dramatic 
changes. If the coupling strength passes over a critical value, a segregation of the decay widths 
occurs. Finally, for large couplings, K rapidly decaying modes and N — K narrow resonances 
result from this interference process @, ||, § . In the following we stick to the statistical approach 
|7|, ||, which is more appropriate for the case of many resonances. Then the internal dynamics 
of the system in consideration is modeled by a matrix ensemble. This may be chosen to describe 
regular and chaotic motion respectively. 

The distortion of regular motion by coupling a closed quantum system to the continuum 
has up to now received little attention theoretically H, |l(], |ll]] and experimentally [12, 13, 14 1. 
In contrast to that, such problems draw considerable attention in classical and semiclassical 
theories (e.g. [15] and references therein). It would be desirable to connect the approaches from 
both fields. The present work is ment to be a first step along that way. 

In [^0|, the regular motion is described by a Poissonian Orthogonal Ensemble (POE). Its 
perturbation due to its coupling to the continuum of decay channels is investigated as a function 
of a parameter a which characterizes the coupling strength. Level repulsion appears at large 
a where the widths have already segregated. There with increasing number of channels the 
correlations in the positions of the resonances approach those, characteristic for the Gaussian 
Orthogonal Ensemble (GOE). At the same time the distribution of the widths experiences a 
considerable broadening in comparison to the Xk distributed widths in the GOE case. This 
broadening is quasi independent of the number of channels. A similar effect has been found in 
the transition-strength dsitribution of a closed system [16, 17] by increasing the chaoticity of 
the system. Note that a broad width distribution implies deviations from the exponential decay 
law g§, ||. 

A physical example of a regular system, coupled to the continuum, is investigated in []TT| , 
13] on experimental and theoretical grounds. The regular motion of the microwaves inside a 



rectangular resonator is disturbed by an attached antenna, which defines the decay channel. 
Theory predicts the effective coupling of the antenna to the resonator to increase with the 
frequency. The experimental data was obtained in the regime of strong external coupling, where 
the resonance widths have already segregated. It should be possible to apply the results obtained 
in the present paper to this experimental setup. 

In the following we compare the results for a regular system described by a POE with those 
of a chaotic system described by a GOE both coupled in the same way to the continuum. Our 
interest focuses on the width distribution and the correlations induced by the coupling of the 
system to the continuum of decay channels. 

In Sec. 2, we describe the statistical model used in our investigations as well as some tech- 
nicalities, concerning the numerical studies performed. A redistribution of the spectroscopic 
properties takes place in the system if its coupling to the continuum of decay channels is suffi- 
ciently strong. We sketch this mechanism (trapping effect) and give the formula for the mean 
width of the long-lived resonances in Sec. 3. The distribution of the widths most characteristic 
of the trapping effect is considered in Sec. 4. In Sec. 5, we present our results for the correlations 
in the positions of the resonances while correlations between widths and positions are discussed 
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in Sec. 6. Some conclusions are drawn in the last section. 



2 Model 

2.1 Projection formalism 

Our analysis proceeds from the statistical model [||, ||]. By using the projection operator tech- 
nique H we couple ./V bound states i = 1,...,N) to K common decay channels 
(\Xc(E)), c = 1, . . . , K) via the two body residual interaction V . We are mainly interested in 
the bound state part of the system and its perturbations by the environment of decay channels. 
To this end, we investigate the statistical properties of the system as a function of the coupling 
strength to the channels. 

The total Hamiltonian of the whole quantum mechanical system consists of three different 
parts: the Hamiltonian H describing the closed system of bound states, a term describing the 
K continua and a third term which specifies the coupling between the two subspaces: 

N K 



#tot = E m H i3 (^-| +J2 J dE \xc(E)) E ( Xc (E)\ (1) 

i,j=l c=l 
K N . 

+EE dE m)vt(E){xc(E)\ + h.c] . 

C=l 1 = 1 

Here the Vf(E) are the coupling matrix elements between the bound states |3>j) and the scatter- 
ing states \xc{E)). They can be understood as the components of a ^-dimensional vector V c . 
Its norm is a measure of the coupling strength to the channel c. Direct reactions are neglected. 
We consider scattering systems with time reversal invariance, so we can restrict ourselves to real 
matrices Hij and Vf(E). 

As shown in || the S-matrix of such a system reads: 

N 

Sab(E) = 5 ab - 2m J2 V?{E)g i3 (E)V b 3 {E) (2) 

where Q is the propagator in the perturbed system H + W xt , 

g = [E - H - W"*] -1 , (3) 
and the generally complex operator W cxt describes the perturbations from outside, 

w , 7= ± i de m^m. (4) 

13 ^ I E + ie-E' 



c=l 



According to (^) , the pole structure of the S-matrix is given by the eigenvalues of the operator 
H + W cxt . We suppose, that the quantities Vf(E') are (almost) independent of the excitation 
energy, which is justified if one considers a finite energy region sufficiently far from the reaction 
thresholds. With this assumption (Q) simplifies 

K 

mr = -^n v ^ (5) 

c=l 
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where we replaced Vf by ^J^-Vi c , Vic being a N x K matrix with on the average unit column 
vectors and rj being the total coupling strength per channel. Then the S-matrix takes the form 

S{E) = l-ir l V + [E-H cB )- 1 V . (6) 

Its complex poles are given by the eigenvalues £ n of the non-hermitian Hamiltonian 

H cS = H - VV + (7) 

with 

£n = E n — -T n . (8) 

Under the influence of the external coupling, the eigenstates of the hermitian matrix H are 
turning into resonances with a finite lifetime. In the case of non-overlapping resonances the real 
part E n and the imaginary part T n in (|8|) give the position and the total decay width of the 
resonance. We will stick to these terms in all cases having in mind the definition as poles of the 
5-matrix. 



2.2 Scattering ensemble 

Following H we represent H eS in the eigenbasis of its hermitean part H. Demanding orthogonal 
invariance of the scattering ensemble, the matrix ensemble is characterized solely by the eigen- 
value distribution of H, and the anti-hermitean part —^VV + is independent on the specific 
choice of the ensemble. In all cases V consists of K random column vectors. For large iV this 
implies the elements of V being independent random Gaussian variables with 



(Vic 



0, 



(V 2 ) = —. 

\ lal 



(9) 



The pole distribution for the .£f off -ensemble may be expressed by 



H eS - Ej + ~Tj 



P{{Ej,Yi}) = c[d N {e k }d NK {V la }l[s(det 

j 

m<n \ n / \ n,a ) 



(10) 



C is a normalization constant. The indices j, k, I, m, n are running from 1 ... AT while a num- 
bers the channels 1 ... K. The integration runs over the whole ensemble parametrised by the 
eigenvalues {e^} of H and the coefficients {Via} of the coupling matrix. By means of the 5— 
functions, the positions and widths of the resonances are introduced as the new variables of the 
distribution. The first part in the last line represents the eigenvalue distribution of H, where a 
is related to the range of the spectrum. The last part in this line represents the distribution of 
the matrix elements of V. 

The parameter (3 > — 1 controls the degree of level repulsion. (3 = 1 refers to the GOE and 
[3 = to the POE case. (3 i— > oo describes a completely rigid spectrum (harmonic oscillator) 
while (3 < refers to a spectrum showing level clustering. Note that the cases (3 = 2 and 
(3 = 4 do not describe the unitary and symplectic ensembles, because we always work with real 
matrices V, demanding orthogonal invariance. 
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In the following we restrict ourselves to the cases (3=1 (GOE) and (3 = (POE). For the 
GOE the level density is described by Wigner's semicircle |2H |, a being its radius. For the POE 
( |lC| ) would produce a Gaussian shaped level density. But in order to compare the correlations 
in both ensembles we find it more convenient to have equal level densities. Therefore we choose 
the {ej} in the POE case as distributed according to Wigner's semicircle law, too. 

Unless stated otherwise, all numerical calculations were done by diagonalising H cS of dimen- 
sion N = 300 and with K = 3 channels being approximately orthogonal to one another, due to 
the centralised distribution of their components. V is a N x K matrix with random Gaussian 
coefficients with zero mean and variance 1/N . The calculations are performed for 50 matrices 
chosen from the GOE or the POE respectively. Then the statistical observables are calculated 
as averages over the central part of the spectrum and over the ensemble (the 50 matrices) si- 
multaneously. In the case of correlation measures (E 2 and Y(q in Sec.|| and Sec.||) the positions 
of the resonances are unfolded to equal mean level spacing d = 1. 

3 Resonance trapping 

3.1 Separation of time scales 

One of the specific properties of open quantum systems is the formation of separate time scales. 
In the strong coupling regime r/ 3> 1, K resonances become very broad by 'trapping' the re- 
maining N — K ones ]|, |2^] . Mathematically the trapping effect is caused by the fact that the 
rank of the hermitean part of H eB is N while that of the anti-hermitean part —^VV + cannot 
be larger than K < N due to its dyadic form (^). Thus we have the following picture: At weak 
coupling the resonances do not overlap: (r) <C D, where D is the mean level distance, and the 
anti-hermitian part of H eS may be treated by first order perturbation theory. Then the widths 
of the resonances are well approximated by the diagonal elements of VV + and have comparable 
magnitude, 

^ = vJ2 V "a- (11) 

a 

In the opposite case of strong coupling, the resonances overlap and the anti-hermitean part 
— ^VV + dominates the behavior of the decaying system. The appropriate basis is now the 
eigenbasis of VV + consisting of K channel vectors (the columns of V) and N — K vectors, 
spanning the kernel of VV + . The channel vectors have the common eigenvalue rj. The widths 
of the remaining states are zero. 

In order to get some information about the magnitude of the widths of the trapped resonances 
it is necessary to go beyond the simple Ansatz sketched above. This is done in the following 
subsection. 

3.2 Mean width of the trapped resonances 

We will calculate the mean width of the trapped resonances in the strong coupling regime 
following H |23|. In addition to what has been done there, we show explicitely the applicability 
to the POE case (in fact it may be applied to any orthogonal invariant ensemble for arbitrary 
[3 as defined in ( JlOj ) ). In order to exploit the orthogonal invariance we need the first and the 
second moments of the distributions of the matrix elements in a generic (random) basis. For 
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the GOE they are pi 



{Hij) = 0, (Hf, 



2N 

,2 



I = J 



' ] a2 . . • ( 12 ) 

AN : * ^ 3 

For the POE we calculate these moments starting from the level distribution, and applying 
a random orthogonal transformation O: 



H = O diag{^}0 T , <Oy>=0, (0%) = jj. 



So every element of H is defined by 



Hij — ^2 OikEkOkj- 



(13) 



(14) 



Due to the centralised distribution of the {Ei} the first moment is zero again. Then we 

calculate the average of the squares Hfj. As the eigenvalues are not correlated, we get 



(Hi) 



Sfc Oik £ k + J2k^i o ik o U EkEi 



is 



3a 2 
4iV 



■IN 



■ i = J 



(15) 



Note that only the diagonal elements differ from the GOE case in their second moment. 

In the following we show that the mean width of the trapped resonances can be calculated 
from the second moments of the matrix elements Hij. To this end we turn to the channel 
representation of H cft . In this basis the channel vectors are the first K canonical vectors and the 
remaining N — K vectors constitute a basis for the kernel of VV + . Then 



H„, 



H 



2V 



(16) 



a, b serve again as channel indices while the greek letters fj,, v are used in connection with 
the long-lived resonances. We diagonalise H in kernel of VV + by the following orthogonal 
transformation 



P 



P, 



D = P H aB P 



( H ab — Y^ab 



V 



'K+l 







z'n I 



(17) 



Note that {e' u } differ from the eigenvalues {ej} of H. D is symmetric and therefore the subma- 
trices denoted by * are transposed to each other, 



N N 

D/j,b = ^k^HklPlb = ^2 Pvp-Hvb- 

k,l=l v=K+l 



(18) 



The Pvp are only correlated with the matrix elements H^ u as they are diagonalising this sub- 
matrix, but not with H^b- Therefore one may average over these quantities independently. 

Now we calculate approximate eigenvalues of D by applying a Jacobi transformation of 



dimension K + l [24]. Such a transformation is trace invariant (this is important, because the 
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total width should be conserved). So we look for the zero of the following determinant (for 
arbitrary K + 1 < fi < N) 



det 



D 



ml 



D 



fj,K 



\ 



D 



Ml 



— H a b + ^5 a b 



Taking into account only the highest orders in rj we arrive at 



0. 



(19) 



4 - £ m + E 



- + O(r?- 2 ) = 0. 



a ° J1 aa T 2 



(20) 



We can read of the approximate position and width of a trapped resonance from the real and 
the imaginary part of For ijwoowe get 



' a 



(21) 



The final step is the calculation of the ensemble mean of D 2 (using (12) and (Il5| 



2 ) 



<£!> 
N ' 



From that it follows that the mean width for the trapped resonances is 



(22) 



(23) 



Note that the diagonal elements of H do not enter in the calculation of (r^), so that ( |23| ) holds 
for both, the GOE and the POE. 

Now one may consider 2^/ (e 2 ) as the 'width' of the level density and then introduce a 
dimensionless parameter k measuring the degree of overlapping: The summed width per channel 
divided by the 'width' of the spectrum, 



For the trapped resonances this overlapping parameter is 

(N-K)(T,) 



Then it follows from (^) that in the limit N ^ oo 



K f 



N 



(24) 



(25) 



(26) 



Assuming ergodicity (as expressed by the randomness of the eigenbasis ( |l3| ) ) the ensemble 
average may be replaced by the spectral average. Then (^) may even be applied to a single 
system. In this case it may serve as a test for the trapping scenario to occur in the system 
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considered: ( |23| ) relates the product of the average width (r^) of the trapped resonances and 
the average width r\ of the broad resonances to the variance of the level density (e 2 ). 

For small k one can still consider the mean width of the N — K smallest resonances in order 
to have Kf well defined globally. Then this quantity will be equal to the mean width of all 
resonances, if k is sufficiently small. Therefore: 

We define the weak coupling regime by k <C 1, so that the upper part of (p7|) is fulfilled, and 
the strong coupling regime by k > 1 so that the lower part becomes valid. In between (critical 
region) the redistribution takes place. 

For illustration, in Fig. 1 the trapping process is portrayed for both the GOE and the POE, 
from the low coupling to the high coupling regime. We plot Kf versus k as defined in (|24|), A 
double-log plot is used in order to demonstrate the proportionality relation (^7|). As the level 
density is the same for the GOE and the POE (cf. Subsec. 2.2), the data points of both the 



GOE and the POE case follow the same line for all coupling strengths within the numerical 
accuracy. The asymmetry on the left side of the plot is due to the finite dimension ./V = 300 in 
our calculations. 

The redistribution between the two scenarios at small and large coupling rj occurs rather 
promptly at k ~ 1. Here the K broad poles appear which will share almost the total sum of 
the widths, whereas the remaining N — K resonances will become more and more long-lived in 
the strong coupling regime. The point k = 1, at which Kf reaches its maximal value, is called 
critical point. Note the peculiar properties of the transmission coefficient [25, |9| and of the width 
distribution [20| at this point. 



4 Width distribution 

In the weak coupling regime the widths are approximately given by the diagonal elements of the 
coupling matrix VV + (11). These are sums of K random Gaussian variables and therefore x\ 



distributed when normalised to unit mean 

= **(*)= rw| yf " le ~ f2/ - (28) 

For the GOE the same distribution holds in the strong coupling region as shown for the one 
channel case in [|[, by calculating the joint probability distribution for the complex eigenvalues 
of H ctt explicitely. For the POE the width distribution becomes much broader. 

The numerical results are shown in Fig. 2. Here a series of width distributions for the GOE 
(left side) and the POE (right side) for increasing coupling strength k (from top to bottom) is 
given. When k > 1, only the long-lived resonances were taken into account. In each diagram, 
we plotted the numerical data (histogram), and the Xk=3 distribution (dashed line). For the 
POE at k = 100 a best fit \\ distribution (dotted line) is given. This distribution is calculated 
by replacing in ( ^8|) A' by a real parameter q and performing a x 2 fit i n order to find the best 
estimate for q. 

For both ensembles the GOE and the POE, the width distributions undergo strong deforma- 
tions in the critical region, but become stationary again in the strong coupling limit. In contrast 



S 



to the GOE case, which returns to its old shape, the distribution for the POE remains much 
broader. It is in good agreement with the best fit \ 2 q distribution with q = 1.7. 

Further investigations of the POE at strong coupling for more channels led to the following 
interesting behavior (cf. Fig. 3): The coupling to continuum leads to a broadening of the width 
distribution. For increasing number of channels K its variance approaches double what it was 
at small coupling (or double what it is in the corresponding GOE case). Whereas the variances 
of the best fit distributions systematically underestimate the factor 2, the numerical variances 
overestimate it. 

In order to understand the broadening of the width distributions due to the coupling to 
continuum, we consider the one channel case in more detail. The S-matrix (|6|) with complex 
argument £ can equivalently be written as [^] 



S{£) 



1-LR 



R 



1 



-V. 



1 + iiT 2 ■ £ -H 

From that it follows that the poles of the S-matrix are given by the zeros of the function 



£ 



Si 



(29) 



(30) 



The Vi are the Gaussian random coefficients of the coupling matrix (11) with variance -k and 
the Si are the real eigenvalues of H. We may consider the real and the imaginary part of f 
separately. As we are interested in the case ij^oowe only keep the highest order of r), making 



use of r ~ rj 1 (f23|). Using the notation £ = E — iT for the poles one arrives at [11] 



E 



E-Ei 



(31) 



Between every two neighboring levels there has to settle down one pole with increasing r) due to 
the structure of (pl|). Choosing an arbitrary one of them and taking into account only its two 
neighboring levels, we find the following formula for its width (a detailed derivation is given in 
App.0): 



-(1-T 2 ), 



V 2 

v 2 



- (32) 
r)(v{ + vi)^ ' n ' v( + v 2 2 

Here s is the distance between the two consecutive levels measured in units of the mean level 
distance do in the center of the spectrum. In our case of a semicircular level density with radius 
a, do = Tfif. In App.[B] we calculate the distribution of the trapped widths normalised by their 
mean (r^) (23) for the GOE and the POE case. This is done by evaluating the following integral: 



p(y 



r/(r,)) = ^-J dsP(s) J d Vl dv 2 e 



7T 2 S 2 1 



4iV vi + vi 



(33) 



Here P(s) is the nearest neighbor distribution for the ensemble in consideration. The result for 
the GOE is: 

1 2y\~ 3 / 2 



and for the POE it is 



pg(v) 



pp(y) 



V^y 



i y 



i + — 

7T 



e^W_ 1>0 (J) 



(34) 
(35) 
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where 

V7T 



(36) 



is the Whittaker function [Eg] and Kq{z) and K\(z) are the modified Bessel functions [g7|. In 
Fig. 4 both distributions are compared to the Porter Thomas distribution Xk=i ( P^l ) an d the 
numerically obtained width distribution for the POE case. Note that we plotted p(ln y) instead 
of p{y) for all distributions in this figure, and we use a logarithmic scale on the abscissa. This 
is done in order to give an overall view and allowing in the same time for the recognition of the 
interesting features discussed in the following. 

Comparing ([34|) to the Porter Thomas distribution, we find that for y <C 1 it agrees exactly 
in the leading power. For large y 3> 1 pc fails: it has a y~ 2 tail which actually makes the 
distribution non normalisable. 

For the POE case we observe the same features when comparing ( |35|) to the width distribution 
obtained numerically. For y — ► oo, W-\<a{z) — > exp(— z/2)/z leading again to a y~ 2 tail being 
inconsistent with the numerical result. On the other hand for y — > 0, pp{y) fits very well to the 
numerical distribution. Here K\{z) — > 1/z and Kq{z) — ► — lnz and therefore 

- (37) 

In the relevant range y ~ 1 both distributions pc{y) and pp(y) show qualitatively how 
the level repulsion parameter /3 ( |To|) affects the width distribution of the trapped resonances. 
Disregarding the tails, one can clearly see from Fig. 4 the broadening in the POE case {(5 = 0) 
in comparisson to the GOE {(5 = 1). 



5 Correlations in the positions of the resonances 

In order to measure the correlations between the positions of the resonances, we apply the 
number variance £ 2 [21]. We investigate the GOE and the POE spectra in the different coupling 
regimes: weak, critical and strong coupling. For the GOE similar investigations have already 
been done in [28]. 

X 2 measures the deviation of the accumulated level density from a straight line. Small S 2 
is a signature for high rigidity of the level sequence. In this case the level spacings are more 
or less equal. On the other hand a completely uncorrelated sequence, as for example the POE 
spectrum has minimal rigidity and therefore maximal S 2 . 

In contrast to |1C] we focus on the few channel case (typically K = 3), because we are mainly 



interested in the differences the GOE and the POE show at large coupling. Note that one result 
of [10] was that the POE spectrum at strong coupling becomes more and more GOE-like with 



increasing number of channels. 

For the numerical calculations we implemented the following technical steps: The spectra 
are unfolded to constant mean level spacing d = 1. This is done by a polynomial fit of the accu- 
mulated level density. Then the edges of the spectrum are skipped, which reduces the number 
of resonances by approximately 25%. Furthermore when k > 1, the K broadest resonances are 
omitted. For the new reduced spectrum {e^} consisting now of N' levels the correlations are 
investigated. 

For an unfolded sequence {e^}, the mean number of levels found in an interval of length I 
is {n{l)) = I due to d = 1. The number variance X 2 (Z) is defined as the variance of n{l) in the 
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ensemble mean: 

T?{l) = {n\l))-{n{l)f = {n\l))-l\ (38) 
For the GOE and the POE (completely uncorrelated sequence) this quantity is known analyti- 



cally [ 21 1 . Instead for our purposes it is sufficient to use the following approximate expressions: 



For the GOE 



X 2 ^)-^ +0.442 (39) 



and for the POE 

Xf>(f) « l. (40) 

In Fig. 5 we show numerical calculations of X 2 for the GOE (a) and the POE (b) spectra. 
This is done for three different values of the overlapping parameter n = 0.1,1, and 10. The 
smallest and the largest value are chosen such that a further decrease or increase of n does not 
change any more the outcome for X 2 significantly. 

X 2 for the GOE (Fig. 5a) shows good agreement with Xq (^) for k = 0.1 and n = 10. Only 
in the critical region we notice a deviation, where X 2 becomes slightly larger for all I > 1. This 
decorrelation is due to the additional 'degree of freedom' the resonance poles encounter in this 
regime [^]. In contrast to the two other cases, where the poles are restricted to a small stripe 
along the real axis, they now have enough space in the direction of the imaginary axis in order 
to avoid close neighbors. This 'width repulsion' cannot be detected by observing the positions 
alone. 

In the POE case (Fig. 5b) all correlations are induced by the coupling to the decay channels. 
Therefore they grow with increasing k. After completion of the redistribution (i.e. when the 
trapped resonances come sufficiently close to the real axis again) the X 2 curve becomes stable. 
This happens at k = 10. Note that the deviations from the original straight line do not vanish 
nor does the curve at large k coincide with that of the GOE. Furthermore it is remarkable, that 
at k < 1 the rigidity increases mainly for large I and catches up at k > 1 for small I. 



6 Correlations connected with the widths 

In [30 1 the number variance is generalized to an intensity variance by weighting each level with 
the intensity of its line as it appears in the cross section. Here we do the weighting with the 
widths of the resonances (in the one channel case both ways are identical). We denote this 
quantity by X 2 (Z). It gives the variance of the summed width in a given interval of length I. 
Neglecting possible correlations with the widths we relate X 2 (Z) to the simple number variance. 
Then we compare the numerical results with this formula. Occurring discrepancies indicate the 
existence of exactly those correlations which had been neglected in the beginning. 
Following |3(J, we define the width-weighted stick spectrum 

Pg( x ) = ~~ e «)> / Pa( x ) dx = N ' ( 41 ) 

with the normalised widths yi = Ti/(T). For k > 1 the K broad resonances are again skipped. 
The summed width in an interval A of length I is 

B{1) = [ Pg (x) dx (42) 
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and consequently the variance of the summed width is 



E 2 (/) = (B 2 (l)) - (B(l)) 2 . (43) 



^2/ 



Now, we relate this quantity to the number variance £ (/) in the following way: 
B 2 ( l ) = ^ViVj 5(x - ei)5(y - ej) dxdy 

v 

= Evil I otherwise +11^ J Six-^-^dxdy. (44) 

i V i^j 
Assuming no position-width correlations, one may perform the averages in ( fl3"l) separately : 

(B(l)) = (y) <n(I)> (45) 
and furthermore assuming no width-width correlations one obtains 

(B 2 (l)) = (y 2 ) (n(0> + (yf E( / S ( x ~ e *)<% " e j) d ^>- ( 46 ) 

It is (n(l)) = I because the mean level distance is one. The corresponding expression without 
width-weighting is 

(n 2 (l)) = (n(l)) + J2([ S(x-ei)S(y-ej)dxdy). (47) 
■ i ■ J A ' 2 

Comparing (|46|) and (|47| ) one arrives at 

( J B 2 (/)) = (y 2 )/ + (y) 2 ((n 2 (0)-/) (48) 
and finally using (|3^) and (||) 



E^(0 = (A^)l + ^(0. (49) 

Here (Ay 2 ) is the variance of the width distribution. In the case that the normalised widths yi 
are Xk distributed, (Ay 2 ) = 2/K. 

From (^9|) we see, that we may clearly separate properties stemming from the correlations in 
the energy spectrum as well as from the width distribution on one hand and correlations between 
energies and widths or between widths at different energies on the other hand. To achieve this 
we proceed |j(J to create from the true width-weighted spectrum (|4l|) a synthetic one 

Pg = ^ZVir(i)5{x - Si) 

where we apply a random permutation it to the indices of the widths. As fl49| ) only uses the 
independence of widths and levels the synthetic spectrum should obey this equation, while the 
true one p g may not if correlations exist. 

Fig. 6 shows the theoretical curves Ep (^) and Eq (^) as well as E 2 for the numerically 
obtained energy spectrum for 3 channels and k = 100. We furthermore show E 2 — (Ay 2 )/ for 
both the true and the synthetic (decorrelated) width-weighted stick spectra under the same 
conditions K = 3 and k = 100. The result for the synthetic spectrum agrees very well with 
that for the energy spectrum as it should be according to (49). The E 2 obtained for the true 
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spectrum on the other hand differs greatly. This is a clear indication that the spectrum in this 
case presents strong correlations, which must be contrasted to the chaotic case where a similar 
analysis indicates the absence of correlations. Note however that the agreement of £ 2 — (Ay 2 )/ 
with Eq for I < 5 is only accidental, what can be concluded from the results for different number 
of channels in Fig. 7. 

The nature of these correlations is such that the spectrum becomes more rigid due to the 
width-weighting (disregarding the term (Ay 2 )/). This may be understood from (32) meaning 
that the width of a trapped resonance is more likely to be large if the distance between its 
neighboring levels is large, too 



Furthermore we expect that for increasing number K of channels the effect of the position - 
width and the width - width correlations vanishes, because the variance of the width distribution 
2/K becomes very small. In this case the sum of the normalised widths becomes more and more 
equal to the number of levels in the interval considered. This is verified in Fig. 7. 



7 Summary 

We considered the distortion of a regular system by its coupling k to the continuum of decay 
channels. We received results for 

«... the width distribution, as it is altered with increasing k: In the weak coupling limit 
for both the GOE and the POE, the widths are x\ distributed when normalised to unit 
mean. Then with increasing k their distributions become broader. For the GOE, the 
widths return to their original distribution in the strong coupling limit. For the POE, the 
width distribution becomes approximately a Xq distribution again, but in contrast to the 
GOE with much larger variance as before (cf. Fig. 2 and Fig. 3). 

«... the correlations between the positions alone: For this analysis we used the number 
variance X 2 . For the GOE we found correlations as they are typical for the level statistics 
of a closed system in both the weak and the strong coupling region. Only in the critical 
region the correlations become weaker. The POE starts without correlations. Then they 
increase steadily with the coupling parameter k (cf. Fig. 5). 

«... the correlations connected with the widths: Here we used the generalized measure S 2 
[j30[l , calculated from a width-weighted stick spectra. There appear two additional types 
of correlations, namely such between the position and the width of each resonance and 
such between the widths of different resonances. In the GOE case they appear only in the 
critical region, whereas in the POE case they increase steadily (cf. Fig. 6 and 7). 

Furthermore we derived an analytic expression for the width distribution of the GOE and 
POE in the one channel case at strong coupling (cf. Fig. 4). It relates the width of a trapped 
state to the distance of the two neighbored levels and explains by this the different width dis- 
tribution of GOE and POE. 



The results of our investigations show the special role of the GOE. Its properties survive the 
distortion of the system by coupling it to the continuum: at large coupling the correlations and 
the width distribution are the same as at low coupling. In contrast to the GOE, the properties 
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of the POE are not restored at strong coupling strength. 

Realistic systems are often in the critical region where correlations in the spectrum are 
induced by the coupling to the continuum. Under these conditions the S-matrix poles are 
difficult to find. Nevertheless they determine the statistical properties of the cross section. We 
will investigate this problem in a forthcoming paper for both the GOE and the POE. 
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A Resonance width in the two level approximation 

Writing ([H]) for two neighboring levels, we get 



+ 



. A 1 _ A 
S ' 2 S 2 



+ 



(s + ±) 2 (a 



(50) 



where A = 62 — S\ is the distance between the levels and E is substituted by E = \{e\ + £2) + s. 
Due to the first equation 

A v\ - v\ A 



2 v\ + v\ 



-T. 



Inserting this result into the second equation of (pOD yields 



A^ 



v i 



+ 



V 2 

' 2 



(l_ r) 2 (1 + r) 



1-T 2 ' 



(51) 



(52) 



Finally measuring A in units of the mean level spacing A = <1qs in the center of the spectrum, 
and r in units of the mean width V = (T u )y according to (p3|), one arrives at ( j32[) 



7T 2 S 2 



y 



4N(v 2 + v 2 ) 

B Width distribution 

The width distribution according to (|32"|) is 



(53) 



p(y) 



N 
2^ 



/J\T / Q Q \ p 2 2 

d»i d» 2 e~T(»l + "2) J i s p( 5 ) { L _ _i£__^_ (1 _ T ^ 



(54) 



Transforming v\ and V2 into spherical coordinates r/yN cos 4> and r/viVsin</> so that r 



sin 



cos 



cos 2(j) we get 



p(y) 



1 

2^ 



rdrdfie 2 / dsP(s)5 



-K Z S Z S1XT 2(f) 
47" 2 



(55) 
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It is enough to integrate <p from to 7r/4 because of the eight-fold symmetry of the integrand. 
Applying in addition the transformation x = sin 2 2<j) 

If —tLf 1 dx f 

In order to go ahead we consider the two cases GOE and POE separately in the following 
subsections. 



y - 



7T 2 S 2 X 

4r 2 



(56) 



B.l GOE case 

Here the nearest neighbor distribution reads 



7T _zu2 



P G (s) = -se-** . (57) 

In order to resolve the S function, we substitute r as a function of y' = ^ ^ x and use / S(y — 
y')f(y')dy = f{y). This leads to 

^=wi' di ^r d " h ~ a< a= ^ i+ ^- <58) 

The integral over the level spacing s gives -l/(2a 2 ) and it remains a last integration, namely 



This may be solved by substituting t = J 



4 r°° t 2 2v 

MV) = -o7T-Tv> / d * . h b = 4 (60) 



vr 2 (l + 6) 2 7o (t 2 + ^r) 2 

and integrating by parts. Then 



B.2 POE case 

Here the nearest neighbor distribution is 

Pp(s)=e- S . (62) 
In contrast to the GOE case we first substitute s instead of r. It follows that 



pp( s ) = — r 2 dre 2 / — e 71 V x 

n z ^J Jo Xy/l — X 

2 f , 2 _ri f°° dze~ az 2r ^ 

= -j— / drr e 2 / , a = — y/y. 

^ \fy J Ji Vz 2 -l vr 
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The last integral represents the modified Bessel function |2 

2 



PP{s) 



drr z e 2 K ( 



This integral can be found in [26] 
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Figure 1: The overlapping parameter Kf of the trapped resonances versus the overlapping pa- 
rameter k of all resonances. The diamonds denote data points for the GOE, the crosses those 
for the POE. The solid lines show the two asymptotics k and k^ 1 . 
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Figure 2: Normalized width distributions for the GOE (left side) and POE (right side) for 
different coupling strengths k = 0.01,0.1,1,10,100 (from top to bottom). Histogram of the 
width distribution obtained numerically (solid line). Xk=3 distribution (dashed line). Best fit 
Xq=i.7 distribution (dotted line) on the bottom left. 
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Figure 3: Relative deviations of the variance of the width distributions (POE, strong coupling) 
for K = 1,3,6,10 channels. The diamonds denote the variances calculated directly from the 
data, whereas the crosses denote 2/q, the variances of the best fit Xq distributions. These values 
were divided by which seems to be the limit value for the variance at K \— ► oo. 




In y 

Figure 4: Distributions of the logarithms of the widths in the strong coupling limit for the 
one channel case. The prediction of the two level approximation for the POE ([35|) (solid bold 
line) is compared to the numerical result (solid line) produced by diagonalising 200 matrices at 
n = 100. The corresponding two level GOE result (|34|) (dashed bold line) is compared to the 
exact distribution, the Porter Thomas curve (dashed line). 
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Figure 5: Number variance S 2 (/) for three different values of the overlapping parameter k = 0.1 
(diamonds), k = 1.0 (crosses) and k = 10. (squares), a) POE and b) GOE. The solid lines show 
the theoretical curves: E 2 G (a) and E£ (b). 
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Figure 6: Variance of the summed width subtracted by (Ay 2 )/ for the POE case. The 

crosses denote this quantity for the original spectrum which contains all correlations. The 
squares correspond to the shuffled spectrum and the diamonds denote the pure number variance 
£ 2 . The solid lines show the theoretical curves Sq and £p. 
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Figure 7: Variance £ 2 (Z) of the summed width subtracted by (Ay 2 )/ but only for the original 
spectrum, containing all correlations. The diamonds denote the case K = 6, the upright crosses 
K = 3, the squares K = 2 and the crosses K = 1. The solid lines show again Sq and Sp. 
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